Tail States in a Superconductor with Magnetic Impurities 
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A field theoretic approacli is developed to investigate the profile and spectrum of sub-gap states 
in a superconductor subject to a weak magnetic impurity potential. Such states are found to be 
associated with inhomogeneous supersymmetry broken instanton configurations of the action. 



While spectral properties of an s-wave superconduc- 
tor are largely unaffected by weak non-magnetic impu- 
rities the pair-breaking effect of magnetic impurities 
leads to the destruction of superconductivity. Remark- 
ably, the suppression of the quasi-particle energy gap is 
more rapid than that of the superconducting order pa- 
rameter, admitting the existence of a narrow 'gapless' su- 
perconducting phase Q. According to the conventional 
(mean- field) theory due to Abrikosov and Gor'kov (AG), 
a gap is maintained up to a critical concentration of mag- 
netic impurities. Yet, being unprotected by the Ander- 
son theorem, it is clear that the gap structure predicted 
by the mean-field theory is untenable and will be de- 
stroyed by rare configurations of the random impurity 
potential. Indeed, since the pioneering work of AG, sev- 
eral authors have explored the nature of the 'sub- 
gap' states. The aim of this work is to investigate quan- 
titatively the spectrum and profile of sub-gap states in 
weakly disordered superconductors. 

In the earliest works 1^-^, attention was focussed on 
the the influence of strong impurities. In the unitar- 
ity limit, it was shown that a single magnetic impurity 
leads to the local suppression of the order parameter 
and creates a bound sub-gap quasi-particle state. For 
a finite impurity concentration, these intra-gap states 
broaden into a band. This mechanism contrasts with 
the AG theory for weak magnetic impurities which 
predicts a gradual suppression of the quasi-particle en- 
ergy gap. Defining C = l/''"s|A|, where |A| represents 
the self-consistent order parameter, and Ts denotes the 
Born scattering time due to magnetic impurities, one 
finds i?gap(Ts) = |A|(1 — ^2/3^3/2^ showiug an onset of 
the gapless region when C = 1 (note Ti = 1 throughout). 
However, even for weak disorder, optimal fluctuations of 
the random potential generate sub-gap states. Extending 
the arguments of Balatsky and Trugman [Q, a fluctuation 
of the random potential which leads to an effective scat- 
tering rate I/t^ in excess of I/ts over a range in excess 
of the coherence length, ^ = (D/|A|)^/'^, induces quasi- 
particle states down to energies Eg^i,p{T'g). These sub- 
gap states are localised, bound to the region or 'droplet' 
where the scattering rate is large, and coupled through 
the proximity effect to the rest of the environment. 

The situation bares comparison with band tail states 
in semi-conductors. Here rare or optimal configurations 
of the random impurity potential generate bound states, 
known as Lifshitz tail states , which extend well below 
the band edge. However, the correspondence is, to some 



extent, superficial: band tail states in semi-conductors 
are typically associated with smoothly varying, nodeless 
wavefunctions. By contrast, the tail states below the 
superconducting gap involve the superposition of states 
around the Fermi level. As such, one expects these states 
to be rapidly oscillating on the scale of the Fermi wave- 
length, but modulated by an envelope which is localised 
on the scale of the coherence length. This difference is 
not incidental. Firstly, unlike the semi-conductor, one 
expects the spectrum of the tails states within the su- 
perconducting gap to be 'universal', independent of the 
nature of the weak impurity distribution, but dependent 
only on the scattering time r^. Secondly, as we will see, 
one can not expect a straightforward extension of exist- 
ing theories p|~pl]| of the Lifshitz tails to describe the 
profile of tail states in the superconductor. 

The aim of this letter is two- fold: firstly we will propose 
a general field theoretic scheme which accommodates the 
AG theory as a mean-field, and whose fiuctuations de- 
termine phase coherence effects in the gapless system. 
Secondly, we will relate the appearance of tail states to 
inhomogeneous supersymmetry breaking instanton con- 
figurations of the action. 

To keep our discussion simple, we will take the 
quenched distribution of magnetic impurities to be 'clas- 
sical' and non-interacting throughout. For practical pur- 
poses, this entails the consideration of structures where 
both the Kondo temperature Q and, more significantly, 
the RKKY induced spin glass temperature is smaller 
than the relevant energy scales of the superconductor. 
The remaining energy scales are arranged in the quasi- 
classical and dirty fimits: ep ^ 1/t ^ (|A|, I/ts) where 
1/t represents the scattering rate associated with non- 
magnetic impurities. The random system we consider is 
specified by the Gor'kov Hamiltonian 
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where Pauli matrices {cr™} and {cf^} operate on par- 
ticle/hole and spin indices respectively. Here S(r) and 
Viy) represent Gaussian 5-correlated magnetic and non- 
magnetic impurity potentials with zero mean and vari- 
ance {JSa{'c)JSji{v'')) g = (67ri^rs)^-^(5'^(r — v'')^^^ and 

(v)V (v'Yiy = (27rz/r)~^5''(r — r') respectively, and 
represents the average density of states (DoS). 

Phase coherence properties of the Hamiltonian (|l|) rely 
on its invariance under symmetry transformations. The 
magnetic impurity potential breaks both time-reversal 
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and spin rotation symmetry leaving only charge conju- 
gation symmetry, 

H = -al"" ® af H'^ af ® al"" . (2) 

When applied to the corresponding Gor'kov Green func- 



tion, G+(e) 



"2 



Gi(-e)crf ® cr^«, the same 
transformation converts an advanced function into a re- 
tarded function implying single-particle interference ef- 
fects as e — > To investigate their influence on 
spectral (and transport) properties of the microscopic 
Hamiltonian it is convenient to cast the problem in the 
form of a functional field integral. Following a standard 
route the Gor'kov Green function can obtained from 
the generating function 



m 



exp 



where ^ represent 16-component superfields with indices 
referencing the particle/hole (ph), spin (sp), charge- 

I 



S[Q]=-^ j dr str 



conjugation (cc), and boson/fermion (bf) spaces, and 
the Pauli matrices {cj''} act in the CC space. The ele- 
ments of the superfields ^' are not independent and ex- 
hibit the symmetry relation ^' = — cr™ ^ a2^j^^ with 
7 = diag(icr5'',-cr'i^)BF- 



As with normal disordered conductors Q , when sub- 
ject to an ensemble average over the random impurity 
distribution, the functional integral over the ^ fields 
can be traded for an integral involving matrix fields Q 
weighted by a non-linear cr-model action. Physically, the 
fields Q, which vary slowly on the scale of the mean- 
free path £, describe the soft modes of density relax- 
ation. When perturbed by the superconducting order 
parameter, the extension of the field integral approach 
is straightforward Taking into account the mag- 

netic impurity potential, the local single-particle DoS can 
be expressed as {iy{r; e))v,s — i'{tr{al^ (^a^^'-' Q{r)))Q/16, 
where (• • ■)q = J^,^-^ DQ (• • .)e-5[0I with 



D{dQ)^ - 4 (ieaf ® - |A|a™) Q 



6Ts 



(3) 



Here D — Vpr/d represents the classical diffusion con- 
stant associated with the non-magnetic impurities. The 
supermatrix fields are subject to the auxiliary symmetry 
condition: Q = O afjQ^-f-^af (g) Although 
the soft mode action is stabilised by the large parameter 
epT ^ 1, the majority of field fluctuations of the ac- 
tion are rendered massive: both the order parameter and 
the magnetic impurity potential lower the symmetry of 
the low-energy theory. To assimilate the effect of these 
terms, and to establish contact with the AG theory, it is 
necessary to explore the saddle-point equation. 

Varying the action with respect to fluctuations of Q, 
subject to the non-linear constraint, one obtains the 
saddle-point or mean-field equation, 



Dd{QdQ) + [Q,te(7 
6ts 



PH 
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With the ansatz: Quf = o""' (T™cosh0 -I- ztT™sinh6', 
where 6 = diag(0i, i0)bf, the saddle-point equation de- 
couples into boson and fermion sectors and takes the form 



dl/^0 + 2i ( cosh 6 



lAI 



sinh6i -Csinh(26l) = 0, (4) 



a result reminiscent of the Usadel equation of quasi- 
classical superconductivity This is no coincidence: 
when subject to an inhomogeneous order parameter, the 
same effective action (||) describes the proximity effect in 
a hybrid normal/superconducting compound |16|. In the 
present context, when combined with the self-consistent 



I 

equation for the order parameter ||l8|, the homogeneous 
form of Eq. coincides with the mean-field equations 
obtained by AG 

The AG solution is not unique: for e 0, the 
saddle-point equation admits an entire manifold of ho- 
mogeneous solutions parameterised by the transforma- 
tions Q = TQuF^"^ where T = IpH (Xi Isp ® t and 
t ~ 7(t^^)-^7^^: soft fluctuations of the fields, which 
are controlled by a non-linear cr-model defined on the 
manifold T e OSp(2|2)/GL(l|l) (symmetry class D in 
the classification of Ref. [Q), control the low-energy, 
long-range properties of the gapless system giving rise 
to unusual localisation and spectral properties. (For a 
comprehensive discussion of the physics of the gapless 
phase, we refer to Refs. . ) 

This completes the formal description of the bulk su- 
perconducting phase. The mean-field solution of the AG 
equation defines the global phase structure of the bulk 
states. Soft fluctuations around the AG mean-field de- 
scribe phase coherence effects due to quantum interfer- 
ence. However, within the present scheme it is not yet 
clear how to accommodate sub-gap states in the gapped 
phase of the AG theory. To identify such states, it is 
necessary to return to the saddle-point equation (^ and 
seek spatially inhomogeneous solutions. We will see that 
such configurations necessarily break supersymmetry. 

To keep our discussion simple, let us focus on a quasi 
one-dimensional geometry. The generalisation to higher 
dimensions follows straightforwardly. To stay firmly 
within the diffusive regime, we require the system size 
L to be much smaller than the localisation length of 
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the normal system ^loc. — vL^D, where Lw denotes the 
cross-section. Furthermore, we focus on the interval near 
to the gapless region (i.e. C ^ 1), where self-consistency 
of the order parameter can be safely neglected. To define 
the inhomogeneous field configurations it is convenient to 
deal not with the saddle-point equation itself, but its 
first integral, (d^/^O)'^ + V{d) = const., where 

V{e) = Ai ^sinh^ - 1^ coshej - Ccosh(2^) , 

represents the effective complex 'potential'. Let us de- 
note by 9ag the values of 9i and i9 at the conventional 
AG saddle-point. From the mean-field DoS, iy{e) = 
V Re cosh ^AG, it is evident that for e < -Egap, Im 6'ag = 
7r/2. The corresponding value of Re 6'ag depends sensi- 
tively on the energy, with Re 6'ag — for e = 0. Taking 
into account the condition that the solution should co- 
incide with ^AG at infinity, one can identify a "bounce" 
solution parameterised by d\ = i7r/2 -I- (/>, with real, 
involving the real potential Vr((/)) = T^(i7r/2 -|- </)) with 
endpoint (j)' such that Vr(0') = Vr(0ag)- 
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FIG. 1. Integration contours for boson-boson and 
fermion-fermion fields in the complex 6 plane. The bounce 
solution for e = (labelled as 'b') is shown schematically. 



As usual, associated with radial fiuctuations around 
the bounce, there exists a zero mode and a negative en- 
ergy mode due to translational invariance of the solu- 
tion. The latter, which necessitates a 7r/2 rotation of the 
corresponding integration contour to follow the line of 
steepest descent (c.f. Ref. [0), has two effects: firstly 
it ensures that the non-perturbative contributions to the 
local DoS are non- vanishing, and secondly, that they are 
positive. Turning to the angular fiuctuations, the sponta- 
neous breaking of supersymmetry is accompanied by the 
appearance of a Grassmann zero mode separated by a gap 
from higher excitations. The zero mode ensures that the 
supersymmetry breaking saddle-point respects the nor- 
malisation condition (2^[0])y_s = and that the local 
DoS is non-vanishing only in the vicinity of the bounce. 




curve) together with the expansion in -Egap — e (dotted curve) 
as determined by Eq. (^. Note that the action vanishes as 
e — > i^Eap. The AG solution for the DoS is shown inset. 



Now integration over the angles 9 is constrained to 
certain contours iQ. Is the bounce solution accessible 
to both? As usual, the contour of integration over the 
boson-boson field 9i includes the entire real axis, while 
for the fermion-fermion field i9 runs along the imaginary 
axis from to ^tt. With a smooth deformation of the in- 
tegration contours, the AG saddle-point is accessible to 
both the angles 9 By contrast, the bounce solution 
and the AG solution can be reached simultaneously by a 
smooth deformation of the integration contour only for 
the boson-boson field 9i (see Fig. |l]). The bounce is as- 
sociated with a spontaneous breaking of supersymmetry 
at the level of the saddle point. 

Having identified the saddle-point field configuration, 
we now turn to the role of fiuctuations. Here we sketch 
the important aspects of the expansion, referring to 
Ref. for a detailed analysis. Field fiuctuations can be 
separated into "radial" and "angular" contributions. The 
former involve fiuctuations of the diagonal elements 0, 
while the latter describe rotations including those Grass- 
mann transformations which mix the bf sector. Both 
classes of fluctuations play a crucial role. 



Taking into account Gaussian fluctuations and zero 
modes, one obtains the non-perturbative, one instanton 
contribution to the sub-gap DoS: 



v,s 



{~i\K\) j z(sinh(/)(a;) — sinh(/)AG)|vo(a;)|^ 
-AttvL^^D\K\sA , (5) 




where S'^ = /^^^ c?0-y/V[t(0AG) — K(0). Here, the factor 

yHTS^jJ^ represents the Jacobian associated with the in- 
troduction of the collective coordinate — is the 
overall factor from the non-zero modes, and (po{x) rep- 
resents the normalised eigenf unction associated with the 
Grassmann zero mode |g5|. Fig. || shows the action 5*0 
for a typical value of C- These results show that, for 
arbitrarily small but finite impurity concentrations, the 
DoS remains finite even at e = 0. For energies e just be- 
low -Egap, an analytical solution can be obtained for the 
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bounce in arbitrary dimension. Generally, one finds the 
exponent Ang{S^/L)'^~'^S^ where 



S^^a^ ^-2/3(1 _^2/3)-(2+<i)/8 



gap 



(6-d)/4 



(6) 



Here g — vDU^^^ denotes the bare conductance, and 
ad is a numerical constant (ai — 8 ^\/24/5). Further- 
more, the optimal solution extends over a length scale 
S.{{Ega,p — e)/\A\)~^/^, set by the coherence length (as 
expected), and diverges as e ^ iJgap (c.f. Ref. ]27|). 

To conclude, let us make some remarks. Firstly, the 
procedure outlined above has a number of close relatives 
in the literature. As well as the investigation of tail states 
in semi-conductors ]Tl| ], a supersymmetric field theory 
was developed by Affleck Q (see also Refs. [^) to in- 
vestigate tail states in the lowest Landau level. There it 
was shown that tail states correspond to supersymmetry 
broken configurations of the ^P-ZieW action (c.f. Ref. [pT[). 
It is also interesting to compare the present scheme with 
the study of 'anomalously locahsed states' Q (see also, 
Ref. There one finds that long-time current relax- 

ation in a disordered wire is also associated with sponta- 
neous breaking of supersymmetry. Finally, we note that 
the Lifshitz argument has been applied on the level of 
the Usadel equation in the study of gap fluctuations due 
to inhomogeneities of the BCS interaction p7| ]. 

Although we have focussed on the question of tail 
states in the superconducting gap, the general scheme 
is more widely applicable. For example, in the present 
system, the transition to bulk superconductivity will 
be preempted by the nucleation of superconducting is- 
lands or droplets within the metallic/insulating phase 
(c.f. Ref. [Q). Similarly, the Stoner transition to a 
bulk itinerant ferromagnet in a disordered system will be 
mediated by the formation of a droplet phase in which 
islands become ferromagnetic |Q. In both cases, we ex- 
pect these 'droplet phases' to be associated with inho- 
mogeneous (replica) symmetry broken saddle-point field 
configurations of the corresponding low-energy action. 
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